Heavy-quark azimuthal momentum correlations as a sensitive probe of thermalization 
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In high-energy nuclear collisions the degree of thermalization at the partonic level is a key issue. 
Due to their large mass, heavy quarks and their possible participation in the collective flow of the 
QCD-medium constitute a powerful probe for thermalization. We present studies with PYTHIA for 
pp collisions at the top LHC energy of ^/s — 14 TeV applying the two-particle transverse momentum 
correlator {Apt,i, Apt, 2) to pairs of heavy-quark hadrons and their semi-leptonic decay products as 
a function of their relative azimuth. Modifications or even the complete absence of initially existing 
, correlations in Pb-|-Pb collisions might indicate thermalization at the partonic level. 
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I. INTRODUCTION 



High-energy nuclear collisions offer the unique opportunity to probe highly excited nuclear matter in the laboratory. 
At sufficiently high temperature and/or energy density hadrons dissolve and quarks and gluons carrying color charge 
become the relevant degrees of freedom. This state of matter is commonly called a Quark Gluon Plasma (QGP). An 
1— I ' essential difference between collisions of leptons or hadrons on the one hand and heavy nuclei on the other hand is the 
^ , development of collectivity in the latter. Collective flow of hadrons, especially the strange hadrons (p and ft, has been 
^ ' experimentally observed at the Relativistic Heavy-Ion Collider (RHIC) (l| suggesting that collectivity significantly 
'■^ ' develops in the early partonic stage, i.e. among quarks and gluons. Presently, the degree of thermalization among 
^ 1 partons remains a key issue. 

r-| ' Heavy-quark hadrons and their observables are of particular interest when addressing thermalization ^3]. In a QGP, 
I— ' chiral symmetry might be partially restored and in that case quark masses should approach the current quark masses. 
^ While these are small compared to the temperature T for up and down quark and comparable to T for the strange 
, quark, they are large compared to T for the heavy quarks (charm and beauty) [3]. Therefore heavy quarks are mostly 
■ created in the early stage of the collision. Further, annihilation of heavy quarks in the QCD-medium is negligible [4j]. 
Thus, heavy quarks probe the entire history of a high-energy nuclear collision. Heavy quarks participate in collective 
motion provided interactions at the partonic level occur frequently. In general, frequent interactions drive a system 
towards local thermal equilibrium. Thus, collective motion of heavy-quark hadrons is a powerful tool when addressing 
. early thermalization of quarks and gluons in high-energy nuclear collisions. 

' In strong interactions heavy quarks are always created together with their anti-quark and are thus correlated. In 
0^ ,. collisions of leptons it has been shown that a hadron containing a charm or beauty quark carries a significant fraction 
of the initial quark momentum ■ Hence, the initial heavy-quark correlations survive the fragmentation process 
^ into hadrons to a large extent and are observable e.g. in the angular distributions of pairs of D- and D-mesons [I,!!]. 

In high-energy collisions of heavy nuclei, frequent interactions among partons (quarks and gluons) of the medium 
and heavy quarks may lead to a significant modification of these initially existing correlations. On the other hand. 



I hadronic interactions at the late stage are insufficient to alter azimuthal correlations of I?D-pairs pJi]. Frequent 
■ " ' interactions distribute and randomize the available (kinetic) energy and finally drive the system, i.e. light quarks 
and gluons, to local thermal equilibrium. To what extent this also happens for heavy quarks is currently a subject 
of discussion An experimental tool to address this question is studied in the current publication. A decrease in 
the strength of heavy quark correlations in high-energy collisions of heavy nuclei as compared to pp collisions would 
indicate early thermalization also of heavy quarks. 

We have employed the Monte Carlo event generator PYTHIA [l2| which implements a fragmentation scheme model 
that reproduces experimentally observed correlations of Z?-mesons at fixed target energies \0\ . To extend calculations 
in PYTHIA beyond leading order, processes contributing at higher orders were calculated using a massless matrix 
element [l^ applying a lower cut-off in the transverse momentum-transfer scale of the underlying hard scattering 
to avoid divergences in the calculated cross section The PYTHIA parameters were subsequently tuned to 

reproduce these next-to- leading order predictions [l3| . For the hadronization of charm and bottom quarks the Lund 
fragmentation scheme was used fl4|. 

Our calculations at leading order (LO) contain flavor creation processes {qq — >■ QQ, gg — ^ QQ) and lead to an 
enhancement at relative azimuth around Acf) « 180° , i.e. the Z?-meson pair is preferentially emitted back-to-back. 
Here, lower-case letters denote light quarks (q) and gluons (g) while upper-case letters denote heavy quarks and anti- 
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quarks {Q,Q). Next-to- leading order (NLO) contributions such as flavor excitation {qQ — qQ,gQ — ^ gQ) and gluon 
splitting (q — > QQ) have a strong dependence on the center-of-mass energy and become dominant at LHC. It was 
noted in [1^ that these higher-order processes do not show pronounced features in azimuth and weaken or destroy 
the azimuthal correlation between charmed hadrons. By choosing momentum cuts, a somewhat enhanced correlation 
was extracted (15i |. 

With these processes expected to dominate at LHC energies, sensitivity to heavy quark thermalization might be lost. 
We investigate whether the pt correlator as a more sensitive measure of correlations would again provide sensitivity to 
thermalization. We introduce the two-particle transverse momentum correlator as a sensitive measure of heavy-quark 
correlations. This method has the following advantages: 

(i) The correlator is sensitive to non-statistical fluctuations, thus carving out any physical correlation. 

(ii) In case of physically uncorrelated candidate-pairs (e.g. background), the extracted value for the correlator 
vanishes, thus providing a reliable baseline. 

(iii) A localization of the observed correlations in transverse momentum space may help to obtain further insight into 
the origin of the observed correlations in relative azimuth. 



II. EMPLOYING THE TWO-PARTICLE TRANSVERSE MOMENTUM CORRELATOR 

We studied two-particle correlations of DD-pa.iis in the (x(pt)i7 x(-Pt)2)-plane, with the cumulative variable x{Pt) 
defined as 

Pt 

x{Pt) = j P{p't)dp'f (1) 



Here, p{pt) is the inclusive pt distribution, normalized to unity, which is obtained from all DD-pairs used in the 
analysis. For the study of two-particle pt correlations, the x(p()-values of DD-paiis {x{Pt)i,x{Pt)2) are filled into 
two-dimensional arrays. In Fig. [U the two-particle correlation function ^^^^^^ is shown for different values of the 
azimuthal separation ^ 4>d ^ <l>j) of the DD-'psii. For small values in A0 we observe a strong positive correlation 




FIG. 1: (Color online) Two-particle correlations as function of (x(pOii x(Pt)2) of 500k DD pairs in different regions of A</> 
integrated over full rapidity for pp collisions at ^/s = 14 TeV as calculated using PYTHIA (v. 6.406). 

at xi ~ X2 ~ 1- This correlation is most pronounced in the high pt-region related to gluon splitting processes. At 
large values of « 180°, a substantial positive correlation in the high p^-region comes from flavor creation processes. 

The occurrence of non-statistical fluctuations of the event- by- event mean transverse momentum Mpt goes along with 
correlations among the transverse momenta of particle pairs [Ts"!. Such correlations were successfully extracted from 
experimental data employing the two-particle transverse momentum correlator [19-23j . For the study of pt correlations 
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between particles of different charge sign like D and D mesons, the correlator is calculated in the following way: 

_ . nev K 

2^fe=l ^^fe ^^k k=l 1=1 j=l 

where pu and ptj are the transverse momentum of the i*'* and j*'* D- and D-meson. Here, the index i runs over all 
particles, while the index j runs over all anti-particles created in a single pp collision. The inclusive mean transverse 
momentum is averaged over all D and £)-mesons, respectively, and is denoted by pt- The total number of DD pairs 
summed over the number of pp collisions n^^ is given by X]fc=i -^k-^k^ with and the number of D and D 
mesons created in a single pp collision. Note that in the present studies we generated one DD pair per pp collision, 
Nj^ = Nj? = 1. In total, we generated 2M (500k) pp coUisions with a DTJ {BB) pair. 

We studied pt correlations and their dependence on azimuthal separation by calculating the correlator in bins of 
A0. In case of independent particle emission, the correlator (Apt^i, Api_2) vanishes. The DD momentum correlator 
(Apt^i , Apt^2) from our simulations of pp collisions at ^/s — 14 TeV is shown in the left panel of Fig. [2] as a function of 
A4>. The error bars reflect the statistical uncertainties from our finite data sample. The correlator has a pronounced 
forward-backward peaked structure. We observe an enhancement at small azimuth from gluon splitting processes, 
while flavor creation of cc-quark pairs leads to an enhanced correlation at backward angles. We have checked that 
flavor excitation processes, involving a larger number of gluons, lead to a rather flat distribution. Also, our studies 
show that the correlations are even stronger at mid-rapidity when compared to full rapidity which can be attributed 
to the harder particle spectrum at mid-rapidity. Integrating the correlator over all azimuth and full rapidity, we get 
(Apt,i, Apt,2) = 0.199 ± 0.006 GeY^/c^. 

In order to account for a possible change in the single particle spectrum when comparing different collision systems 
or energies, the normalized dynamical fluctuation 

T,pt = sgn((Apt,i, Apt,2)) = (3) 

has been introduced as a dimensionless measure [20j . 

Our result for DD mesons as given above amounts to Ept w 28% with p[^ — 1.58 GeV/c. This represents a large 
value implying a strong correlation when compared to e.g. « 1% observed for unidentified charged particles 
in central collisions at SPS and RHIC |20l - [23| . To mimic combinatorial background which is always present in 
the experiment, we applied the correlator to D- and D- mesons from different pp collisions, which are physically 
uncorrelated. This results in a value of {Apt.i, Apt. 2) consistent with zero (see Fig. [2]). Therefore the correlator allows 
for a clear distinction between the case were correlations are present (different from zero) or absent (equal to zero) 
going beyond the method described in [10|, il5|] . 
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FIG. 2; (Color online) Distribution of the momentum correlator { Apt, 1, Apt, 2) of 500k DD pairs (left panel) and 500k BB 
pairs (right panel) as a function of relative azimuth A(j> at mid-rapidity (circles), for full rapidity (squares) and for background 
using the mixed event method (triangles) for pp collisions at y's — 14 TeV as calculated using PYTHIA (v. 6.406). The lines 
connect the points. 
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FIG. 3: (Color online) Distribution of the momentum correlator (Apt, i, Apt, 2) of 200k e^e~-pairs from DD decays and De 
correlations (left panel) and of 50k e^e~ -pairs from BB decays (right panel) as a function of relative azimuth Acj) for full 
rapidity for pp collisions at ^/s = 14 TeV as calculated using PYTHIA (v. 6.406). The lines connect the points. 

The BB momentum correlator is shown in the right panel of Fig. [2] and has a structure similar to the one for 
UD-pairs. Integrating the correlator over all azimuth and full rapidity, we get (Apt. 1, Apt, 2) = 2.73 ± 0.05 GeV^/c^ 
which corresponds to the normalized fluctuation Spt ~ 31% with = 5.27 GeV/c. This demonstrates that by 
applying the momentum correlator to pairs of heavy quarks in pp collisions at LHC energies, strong correlations are 
predicted which should be experimentally observable. When only considering DD production yields as a function of 
relative azimuth, a weaker dependency is predicted (isj . 

At LHC energies, the production of D-mesons is dominated by gluon splitting and flavor excitation processes while 
the contribution from flavor creation is about 10% at low momentum and increases up to 20% at larger momentum p^l . 
On the other hand, the production of i?-mesons is dominated by flavor creation and flavor excitation with a small 
contribution below 10% from gluon splitting and an overall weak dependence on transverse momentum j2^ . 

As shown above, the initial correlations of cc-quark pairs survive the fragmentation process into hadrons to a 
large extent. However experimentally, full kinematic reconstruction of I?-mesons from topological decays suffer from 
small branching ratios and rather small reconstructing efficiencies resulting in low statistics, especially when pairs of 
D-mesons are considered where these factors enter quadratically. In minimum bias pp collisions at y/s = 14 TeV, 
roughly 28 (2) out of 1000 collisions create a charmed (beauty) meson such as D°,D+ or Z?+ {B^,B^ or B^) at 
mid-rapidity |y| < 1. The branching ratio in the golden channel D^ — )■ K~ + 7r+ amounts to 3.83% with additional 
penalty factors due to the detector acceptance for the decay particles and topological reconstruction of the secondary 
decay vertex. Overall we estimate the number of fully reconstructed D^ — D^ pairs in 10^ minimum bias pp collisions, 
which is equivalent to one year of ALICE data taking, to be in the order of 10. This is obviously too low to study pt 
correlations. 

As an alternative, we considered electrons (positrons) from semi-leptonic decays of charm and beauty hadrons with 
an average branching ratio to electrons of 10% and 11%, respectively. 

The left panel of Fig. |3]shows the momentum correlator versus relative azimuth for D — e-pairs (triangles), with the 
electron stemming from the semi-leptonic decay of one of the D-mesons and for — e~-pairs (squares) where both 
_D-mesons decayed into an electron. The right panel shows the correlator for e"'"e~-pairs from semi-leptonic decays of 
-Bi?-mesons pairs. The correlations at small values of Acf) do not survive the semi-leptonic decay while at backward 
angles around Ac/) w 180°, we still observe a strong correlation. We checked that this is due to the decay kinematics 
with gluon splitting processes dominating at forward angles resulting in a softer distribution of the heavy-quark 
hadron. At backward angles flavor creation processes dominate leading to a significantly harder spectrum. 

Integrating the correlator over all azimuth and full rapidity, we extract (Apt, 1, Apt, 2) — 0.007 ± 0.001 GeV^/c^ 
(0.10 ± 0.01 GeV^/c^) for e+-e~-pairs from charm (bottom) decays corresponding to the normalized dynamical 
fluctuation Ept w 17%(18%) with p[ = 0.50 GeV/c (1.73 GeV/c). This clearly indicates that the initial correlations 
among a heavy quark and its corresponding anti-quark even survive semi-leptonic decays into electrons (positrons) to 
a lar ge e xtent. With the predicted charm production in full rapidity of 0.16 per minimum bias pp collision at -^/s = 14 
TeV [135, estimate the number of electron-positron pairs in the ALICE central barrel {\y\ < 0.9, pt > 0.2 GeV/c) 
from heavy-quark decays within one nominal year of ALICE running to be more than 100k. Thus, an experimental 
observation of heavy-quark momentum correlation at the LHC should be possible. 

The heavy-quark transverse momentum depends on the relative contributions from different QCD-process. Fur- 
thermore, experimentally it has been observed that the average transverse momentum is monotonically rising with 
the charged-particle multiplicity in pp collisions [25l [26j . We studied the normalized dynamical fluctuation Ept of 
DD-paiis in several multiplicity regions in pp collisions at the top LHC energy as shown in Fig.|4l Higher multiplicity 
collisions result in stronger correlations due to the increase of the mean transverse momentum from p7 — 1-53 GeV/c 
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FIG. 4: (Color online) The normalized dynamical fluctuation Ept of _DD-pairs as function of A(j) for different multiplicity classes 
at full rapidity from pp collisions at ^/s = 14 TeV as calculated using PYTHIA (v. 6.406). For each multiplicity class, 10k 
-DD-pairs were used. The lines connect the points. 



at multiplicities {Nch) ~ 25 to pT = 3.46 GeV/c at (Nch) ~ 42. In addition, the contribution from flavor creation to 
the production of D-mesons increases up to 18% in high multiplicity events compared with 12% in low multiplicity 
events leading to enhanced correlations at large relative azimuth A(j3. Thus, the highest multiplicities in pp collisions 
at LHC energies might be a good case to experimentally establish the existence of these correlations for heavy quarks. 

The results on the heavy-quark correlator discussed above are our prediction for pp collisions at the top LHC energy 
and serve as a baseline for the case that no thermalization sets in as is expected for such a small collision system. 
Further, our calculations show that higher multiplicity pp collisions result in stronger correlations. 

Finally, we consider the relative pseudo-rapidity Arj of DD pairs and show the correlator as a function of 
and Arj (see Fig. [S]). Gluon splitting processes lead to correlations at small values of A(p and Arj. On the other 
hand, flavor creation results in correlations at large values of Acj) extending over a large range in A77. In Pb-|-Pb 
collisions, the development of strong transverse flow would lead to a broadening of the away-side momentum correlation 
and an enhancement at small relative azimuth of DD pairs j27l . [2^ . To experimentally disentangle these different 
contributions, an analysis in ranges of pseudo-rapidity, e.g. jAryj < 0.5 to study gluon splitting processes and effects of 
collective flow versus | Aryj > 0.5 where flavor creation dominates, might help. Further, a comparison of experimental 
heavy-quark correlations from Pb-|-Pb collisions to results from microscopic transport calculations would provide an 
independent way to extract effective heavy-quark scattering rates in the QCD-medium [2^ . 



III. CONCLUSIONS AND OUTLOOK 



In summary, we have presented a sensitive method to see azimuthal correlations of heavy-quarks in pp collisions at 
LHC energies. We applied the momentum correlator to pairs of heavy-quark hadrons and their scmi-leptonic decay 
products as a precise and normalized measure. At LHC energies, the production of charm quarks is expected to be 
dominated by gluon splitting processes resulting in forward correlations with an increasing contribution of backward- 
peaked pair creation at larger momentum. A stronger correlation is expected in high-multiplicity pp collisions. A 
modification or disappearance of these momentum correlations in Pb-|-Pb collisions as compared to pp collisions will 
be explored as a sensitive probe of thermalization. 
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FIG. 5: (Color online) The momentum correlator of 2 million DD-pairs as a function of Ac/) and Ar; from pp collisions at y^s 
= 14 TeV as calculated using PYTHIA (v. 6.406). 
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